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ABSTRACT 

We quantify the consequences of intergalactic dust produced by the first Type II 
supernovae in the universe. The fraction of gas converted into stars is calibrated based 



on the observed C/H ratio in the intergalactic medium at z = 3, assuming a Scalo 
mass function for the stars. The associated dust absorbs starlight energy and emits it 
at longer wavelengths. For a uniform mix of metals and dust with the intergalactic gas, 
we find that the dust distorts the microwave background spectrum by a y-parameter 
in the range (0.02-2) x 10 _5 (Msn/0.3Mq), where Msn is the average mass of dust 
£D • produced per supernova. The opacity of intergalactic dust to infrared sources at 

redshifts z ^ 10 is significant, T^ ust = (0.1-1) x (Msn/0.3Mq), and could be detected 
■^■j- ■ with the Next Generation Space Telescope. Although dust suppresses the Lya emission 

from early sources, the redshifts of star clusters at z = 10-35 can be easily inferred 
| from the Lyman-limit break in their infrared spectrum between 1-3.5/iin. 

q ■ Subject headings: cosmology: theory - cosmic background radiation 
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s_j ■ 1. Introduction 

a, 

Recent spectroscopic observations of the Lya forest at z J> 3 show evidence for a 
metallicity ~ 1%Zq in absorbers with HI column densities as low as 10 15 cm -2 (Cowie 1996; 
Songaila & Cowie 1996; Tytler et al. 1995). Numerical simulations identify such absorbers with 
mildly overdense regions in the intergalactic medium, out of which nonlinear objects such as 
galaxies condense (Hellsten et al. 1997, and references therein). Indeed, damped Lya absorbers, 
which are thought to be the progenitors of present-day galaxies, show similar metalicities during 
the early phase of their formation at 3.5 ^ z <^ 4.5 (Lu, Sargent, & Barlow 1996). The universality 
of this mean metal abundance in absorbers spanning a range of six orders of magnitude in HI 
column density, indicates that an early phase of metal enrichment occurred throughout the 
universe at z ^ 5, before galaxies were assembled. 



The formation of an early population of star clusters at redshifts ~ 10-30 is a natural 
consequence of hierarchical structure formation in Cold Dark Matter (CDM) cosmologies (see, e.g. 
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Haiman & Loeb 1997, hereafter HL97; Ostriker & Gnedin 1996, and references therein). The first 
star clusters that form in high abundance have baryonic masses ~1O 8 M0 and virial temperatures 
~10 4 K, which allow their gas to cool via atomic transitions. Less massive objects with lower virial 
temperatures are unable to cool and fragment into stars due to the rapid photo-dissociation of their 
molecular hydrogen (Haiman, Rees, & Loeb 1996). Since the potential wells of the first clusters are 
relatively shallow (~10 km s _1 ), supernova-driven winds might have expelled the metal-rich gas 
out of these systems and mixed it with the intergalactic medium. Incomplete mixing could have 
led to the observed order-of-magnitude scatter in the C/H ratio along lines-of-sight to different 
quasars (Rauch, Haehnelt, & Steinmetz 1996; Hellsten et al. 1997). It is an interesting coincidence 
that the supernova energy output associated with a metal enrichment of ~ 1%Zq corresponds to 
~ 10 eV per hydrogen atom, which is just above the binding energy of these early star clusters. 
Supernova feedback in these objects could have therefore determined the average metallicity 
observed in the Lya forest. Direct observations of these supernovae might be feasible in the future 
(Miralda-Escude & Rees 1997). 

The measured C/H ratio can be used to calibrate the net fraction of gas which is converted 
into stars by a redshift z ~ 3. In HL97, we have used this fraction and an extension of the 
Press-Schechter formalism to calculate the early star formation history in a variety of CDM 
cosmologies. For a large range of models we have found that the universe is reionized by a 
redshift z = 10-20, and that the resulting optical depth to electron scattering, ~ 0.1, is detectable 
with future microwave anisotropy experiments. In addition, deep imaging with future infrared 
telescopes, such as the Space Infrared Telescope Facility (SIRTF) or the Next Generation Space 
Telescope (NGST), would be able to detect pre-galactic star clusters. In particular, NGST should 
find ^ 10 3 star clusters at z > 10 within its field of view of 4' x 4', given its planned detection 
threshold of 1 nJy in the wavelength range of 1-3.5 fim (Mather &; Stockman 1996). 

The early epoch of star formation and metal enrichment is inevitably accompanied by the 
formation of dust in supernova shells. This dust has two important observational signatures. 
First, the absorption of starlight energy and its re-emission at long wavelengths distorts the 
spectrum of the cosmic microwave background (CMB) radiation (Wright et al. 1994; Bond, Carr 
& Hogan 1991; Adams et al. 1989; Wright 1981). Second, the opacity of the intergalactic medium 
to infrared sources at redshifts z ;> 10 could be significant. For these redshifts, infrared in the 
observer frame corresponds to UV in the source frame-a spectral regime in which dust absorption 
peaks.f] Dust obscuration must therefore be considered when predicting the performance of future 
infrared telescopes such as NGST. 

In this Letter, we quantify the magnitude of the above signatures of intergalactic dust. In 
§ 2 we describe our method of calculation and in § 3 we compute the distortion of the microwave 



The effect of intergalactic dust at much lower redshifts was discussed in the context of the reddening (Cheney 
& Rowan-Robinson 1981; Wright 1981, 1990) and the number count (Ostriker & Heisler 1984; Ostriker, Vogeley, & 
York 1990) of quasars. 
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background spectrum and compare it to the current COBE limit. In § 4 we evaluate the opacity of 
the intergalactic medium to infrared sources. Finally, § 5 summarizes the implications of this work 
for future microwave-background and infrared observations. Throughout the paper, we adopt a 
Hubble constant Hq = 50 km sec" 1 Mpc" 1 and a cosmological density parameter Q = 1. 



2. Method of Calculation 

Our first goal is to calculate the distortion of the CMB spectrum due to the production 
of radiation by stars and its processing through dust. The CMB spectrum is initially a pure 
black-body with a temperature Tcmb = 2.728(1 + z) K (Fixsen et al. 1996). We assume that 
stars are born with the Galactic initial mass function (IMF) as parameterized by Scalo (1986), 
and compute their time-dependent spectra using standard spectral atlases (Kurucz 1993) and 
evolutionary tracks (Schaller et al. 1992). Our stellar modeling is described in detail in HL97. 
In computing the amount of dust produced by the stars, we assume that each type II supernova 
yields O.3M of dust which gets uniformly distributed in the intergalactic medium (see discussion 
in § 3), and that the dust absorption follows the wavelength-dependent opacity of Galactic dust. 

To describe the effect of dust on the CMB we follow Wright (1981), and define the comoving 
number density of photons with a comoving frequency v as 

N "( z ) = hc{l + zf Jv{1+z){ ■ Z), (1) 

where J v (z) is the specific intensity in erg cm" 2 s _1 Hz" 1 sr" 1 . The evolution of N v {z) is 
determined by the radiative transfer equation: 

dN v (z) cdt 



dz dz 



j v {z) - a u{1+z) (z)N u (z) , (2) 



where in an U = 1 universe, (cdt/dz) = (c/Hq)(1 + z) 5//2 . Here, j v {z) is the comoving emission 
coefficient due to stars and dust, 
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UZ) = h^TW jl{1+z){Z) + * ie X p M l + ^)ABT dust ]-l) a "(^) (3) 
and a v (z) is the dust absorption coefficient, 

0L V {Z) = p&ust{z)K u , (4) 

for a mass density, pdust) in dust. We assume that the dust is in thermal equilibrium with the 
background radiation field and equate the power it absorbs in the UV to the power it emits in the 
infrared. The dust temperature Td us t as a function of redshift is then derived from the implicit 
equation, 

iVK "^ \ NAA - I exp[Ml + - 1 } = °- (5) 



Jo 
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We adopt the dust opacity k v based on the mean Galactic extinction law A(X)/A(J) of Mathis 
(1990), 



8.3 x 10 3 



-4(A) 
A(J) 



> : u)si „,.,., ) cm2g - (6 ) 



0.12 + (A(/im) -0.4) 2 

where A = cjv is the wavelength; the numerical coefficient converts the extinction per hydrogen 
atom, with the assumed gas-to-dust ratio of 100, to extinction per gram of dust; and the last 
term is the albedo term given by Wright (1981) for converting the extinction (absorption plus 
scattering) to pure absorption. 

Given a specific star-formation history, stellar spectra, and dust-yields, one can integrate 
numerically equations (||H([|) and get the dust temperature and the net spectrum of the radiation 
background at z=0. As in HL97, we express the star-formation history in terms of the mass 
fraction of all baryons which are assembled into collapsed objects, F co \\(z), and the fraction of 
these collapsed baryons which get incorporated into stars, /star- The density parameter of stars is 
then, f2 sta r(-z) = fsta,rF co \\(z)^, where fib is the total baryonic density parameter (= 0.05 in our 
standard model). We adopt the values of both F co \\(z) and / star from HL97; F co u(z) is calculated 
in a standard CDM cosmology based on an extension of the Press-Schechter formalism that takes 
into account pressure and photo-dissociation of H2, and / s t a r ~ 4% is found from the requirement 
that the C/H ratio be 1% of the solar value at z = 3. The stellar emission coefficient is found at 
each redshift through a convolution of the time-dependent emission from all stars born prior to 
that redshift 

&{z) = fstA Pc0 (l + z) 3 J^dz'^e*[u(l + z'),t Z:Z/ ], (7) 

where p C Q is the current critical density of the universe, t ZtZ i is the time interval between redshifts 
z' and z, and e*(v, t) is the composite physical emissivity (in erg s _1 Hz -1 Mq 1 ) of a cluster of 
stars with the Scalo IMF at time t after its initial starburst (see HL97 for details). We have 
assumed that all stellar photons escape their parent clouds and are ejected into the IGM. In reality, 
a fraction of the ionizing photons will get converted into line photons with energies E < 13.6eV by 
case B recombinations inside the clouds, but this will only have a minor effect on the total amount 
of energy absorbed by dust. Free- free emission by the ionized gas in the cluster can be neglected 
relative to the stellar emission. The mass density of dust at each redshift is given by 

Pdust(z) = /dust/starfibPco(l + zf J dz fdep(t Z ,z'), (8) 

where /dust * s the mass fraction of stars which gets converted into dust, and < /dep(^) < 1 is the 
fraction of the total mass in dust which is deposited after a time t following the starburst (i.e. the 
mass fraction of the > 8Mq stars which have completed their main sequence lifetime). Assuming 
that each supernova in a Scalo mass function produces O.3M of dust leads to /dust = 0.001. The 
uncertainties in this value will be discussed in the next section. 
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3. Distortion of the Microwave Background Spectrum due to Dust Emission 

To parameterize the deviation of the spectrum at z=0 from a Tcmb=2.728 K black-body, we 
compute the Compton y-parameter y c , defined as 



1 

Vc = 7 



J dvN v l 
I duN 0iV 



(9) 



where the integrals are evaluated over the FIRAS frequency range of 60-600 GHz (Fixsen et al. 
1996), and 

c 6 exp(/iz//fcB-r C MB) - 1 

The y-parameter we derive is subject to a number of uncertainties associated with the 
formation of dust. Each of these uncertainties can be parameterized by a separate factor that 
multiplies our result for y c . The first uncertainty involves the average fraction of the total stellar 
mass which gets converted into dust, /dust- The value /dust = 0.001 adopted in § 2 is equivalent 
to assuming that each star with mass above 8M© explodes as a supernova and produces O.3M of 
dust. This mass is consistent with the values inferred from observations of the depletion of the 
Si I line in SN 1987A (Lucy et al. 1991), and from the infrared echo of SN 1979C and SN 1980K 
(Dwek 1983). Theoretical estimates with 100% condensation efficiency predict values as high as 
1M of dust per supernova (Dwek 1988; Weaver & Woosley 1980). If the dust yield of supernovae 
is supplemented by a significant amount of dust produced in stars, the value of /dust could increase 
even further. On the other hand, a substantial fraction of the dust may be destroyed by sputtering 
inside the supernova envelope, or by shock waves from other supernovae (Dwek 1997); the 
y-parameter scales linearly with the survival fraction of dust, / surv . The y-parameter would also 
be affected by a spatial correlation between the distributions of dust and stellar sources. If dust is 
correlated with star forming regions, the radiation flux probed by the dust could be enhanced by a 
factor /flux relative to the average flux in equation (0). For small values of /flux, Vc scales linearly 
with /flux- The final uncertainty involves the mixing efficiency of carbon with the intergalactic gas. 
The star-formation efficiency, /stan which is required in order to get the observed metallicity of 
the Lya forest is inversely proportional to this mixing efficiency. From equations (|2|), (Q) and (|7]) 
we get y c oc / s 2 tar , since the radiation emitted by dust is proportional to its opacity and the stellar 
radiation background - both of which are independently proportional to /star- All of the above 
uncertainties can be combined through the definition of an overall efficiency parameter r] y : 

r p ( /dust \ / /star \ f 1 1 *\ 

Vy = /surv/flux {-^) [^) K 11 ) 

so that y c oc r/ y , and rj y = 1 in our standard model. 

In addition to the uncertainties about the dust and the star-formation efficiency, we note 
that if an early generation of low-luminosity quasars of the type discussed by Eisenstein & Loeb 
(1995) exists, these quasars could have produced UV radiation exceeding the stellar output. The 
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UV light from these quasars would be reprocessed by dust in much the same way as the stellar 
radiation is, and could substantially increase the spectral distortion derived from the stars alone. 
The possible effects of early quasars are ignored in this Letter and will be quantified in a future 
paper. 

Figure 1 shows the full spectrum of the background radiation for our standard model, whose 
parameters are summarized in Table 1. The solid line depicts the total spectrum, which is a sum 
over the contributions from the original CMB, the stars, and the dust, shown separately by the 
dotted and dashed lines. The net contribution from dust is negative (absorption) for A & 0.015 
cm. The dust emission at longer wavelengths peaks around A « 0.1 cm, within the COBE range. 
The y-parameter derived from this spectrum is 0.46 x 10~ 5 , a factor of 3 below the upper limit of 
1.5 x 10" 5 established by COBE (Fixsen et al. 1996). 

For reference, the bottom panel of Figure 1 shows the star-formation history, i.e. the fraction 
of baryons confined in stars as a function of redshift, for our standard model. The same panel 
also shows the histories of the dust temperature and the y-parameter. The dust temperature 
starts to deviate noticeably from the CMB temperature (dashed line) around z 10, when the 
fraction of baryons in stars fi s t ar /^b = /star-Fcoii ~ 1%- Note that we ignore any contribution to 
the distortion from redshifts z <3. Our calculation is inadequate at low redshifts, since it ignores 
the enhanced production of metals and dust inside galaxies and the contribution of quasars to the 
background radiation. The spectral distortion we obtain should therefore be regarded as a lower 
bound; inclusion of the contribution from lower redshifts could only increase its magnitude. 

Table 1 shows the range of values obtained for the y-parameter as one deviates from the 
input parameters assumed in our standard model. In particular, if the normalization of the 
standard CDM power-spectrum a 8h -i is increased, then star-formation starts earlier and spreads 
over a wider range of redshifts. Although the optical depth to dust is increased by 70%, the 
corresponding change in the y-parameter is only 10%, since the average redshift delay between 
the production of radiation and dust is increased. Tilting the power spectrum to have a lower 
power-law index n has the opposite effect; stars are produced later and over a narrower range of 
redshifts. In this case, the optical depth is reduced by 50% but the y-parameter decreases only 
by 10%. Since we normalize the star-formation efficiency by fixing the total carbon yield of the 
stars, any change to the slope of the IMF has little effect on the net production of dust. Adding a 
constant —0.7 < (5 < 1 to the power-law index of the IMF changes the y-parameter in the range 
0-17 < (y c /10~ 5 r/ y ) < 0.54. When the photo-dissociation of H2 is ignored (e.g., due to the renewed 
formation of H2 on dust grains), the threshold for fragmentation and star formation in virialized 
objects is lowered from virial temperatures T v - ir ^ 10 4 K (where atomic line cooling is effective) to 
T v ir ^ 10 2 K (where H2 cooling operates). In this case, star-formation starts earlier and leads to 
a slight increase in the y-parameter. The input parameter that affects y c the most is the baryon 
density, y c cx fij^; high values of fib (such as Ob ^ 0.09, if all other parameters are held fixed) 
violate the existing COBE constraint. 
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4. Optical Depth to Infrared Sources 

Next, we evaluate the optical depth to absorption and scattering by dust along the line of 
sight to a source at a redshift z s > 3, 

c f Zs dz 

T dw t{v) = —J 3 {1 + z)5/2 a»(i+z)(z)- (12) 

The absorption coefficient QWi+s) (z) is given in equations (|j) and Since we are now interested 
in dust obscuration along a specific line of sight, both absorption and scattering must be included, 
and so we omit the last term in equation (^). The uncertainties discussed in § 2 concerning the 
formation and survival of dust, as well as the fraction of the baryons converted into stars, affect 
the optical depth through the parameter r] T , 

p p ( /dust \ / /star\ /'IQA 

Vr = /surv/flux ) ) > W 

so that Tdust obtains the value predicted by our standard model times rj T . 

Figure 2 shows the resulting optical depth in our standard model as a function of observed 
wavelength A, for three different source redshifts. The optical depth peaks at A ~ 1/xm, where it 
reaches a value Td us t = 0.4. This high value could suppress the Lya emission peak in the spectra 
of early star clusters (see Chariot & Fall 1991); attempts to detect this emission line from galaxies 
at z ~ 2-4 proved difficult in the past (Thompson, Djorgovski & Trauger 1995; Lowenthal et al. 
1997). Since the total amount of dust is the same in all models with the same IMF, all other 
variations in the peak optical depths given in Table 1 arise from differences in the production 
history of dust. In particular, the optical depth is high when stars form early (e.g., as a result of 
a high Cg/j-i or the lack of H2 destruction), and low when stars form late (low n). We note that 
the optical depths we derive are below the values that could cause a measurable reddening in the 
spectra of quasars at z <J 5 (Rowan-Robinson 1995). 

To illustrate the potential effect of dust on future observations with NGST, we show the 
composite spectrum of two star-clusters, located at z s =10 and 20, and containing 10 8 M Q and 
4.6 x 10 8 M in stars, respectively (bottom panel). The stars are assumed to be distributed on 
the main sequence according to the Scalo IMF. In both cases the top curve shows the original 
unprocessed spectrum, and the bottom curve shows the obscured spectrum due to the intervening 
dust in the standard model (with rj y = 1). These results imply that the performance of NGST 
near its planned detection threshold of 1 nJy in the wavelength range 1-3.5 /im, could be affected 
by intergalactic dust. 

Another important conclusion that can be drawn from Figure 2 is that despite the obscuration 
by dust, the Lyman-limit break could easily be identified in the emission spectrum of star clusters 
at high-redshifts. Additional absorption by neutral hydrogen in the IGM would only make the 
break appear sharper. For NGST observations between l-3.5//m, the detection of this break can 
identify source redshifts in the range z s =10-35. The feasibility of this UV dropout technique was 
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recently demonstrated through the successful photometric identification of galaxy redshifts in the 
interval 2.5 £ z s £ 3.5 (Steidel et al. 1996; Madau et al. 1996). 

5. Summary 

We have quantified the imprint of intergalactic dust due to the first stars on the y-parameter 
of the CMB spectrum (in the COBE wavelength range), y c , and on the opacity to absorption and 
scattering of radiation from sources at high redshifts, Tdust^s)- We normalized the star-formation 
efficiency so as to get the observed ~1% solar C/H ratio in Lya absorbers at z=3. With this 
normalization, we find y c ~ 10 -5 and id us t(oo) ~ 1 for plausible choices of the cosmological 
parameters and stellar properties at high redshifts (see Table 1). Uncertainties due to the mixing 
efficiency of metals with the IGM, the stellar IMF, the destruction of dust, the spatial correlations 
between dust and stellar sources, or the cosmological parameters, could change these estimates by 
an order of magnitude. If early quasars exist, their additional UV radiation would enhance this 
spectral distortion. 

The current COBE limit, y c < 1.5 x 10~ 5 , already places interesting constraints on the dust 
formation history of the universe at z ;> 10. Complementary microwave anisotropy limits on 
the scale of ~ 1" would constrain the degree of patchiness in the dust distribution. But most 
importantly, the future detection of source reddening at z s ^ 10 by NGST, would be instrumental 
in establishing the properties of the first stars in the universe. Despite the significant dust opacity, 
Figure 2 shows that infrared photometry by NGST could identify the redshift of sources at 
10 & z s ^ 35, based on their pronounced Lyman-limit break at observed wavelengths of 1-3.5/um. 

We thank Dimitar Sasselov and Eli Dwek for useful discussions. This work was supported in 
part by the NASA ATP grant NAG5-3085 and the Harvard Milton fund. 
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Table 1: Parameter values in our standard model and its variants, and their effect on the expected 
values of the y-parameter and peak opacity of the intergalactic dust to a source at z = oo. The 
efficiency parameters r] y and r\ T are defined in equations (11) and (13). Models with a predicted 
y-parameter in excess of the COBE upper limit, y c < 1.5 x 10 -5 , are ruled out. The optical depth 
of intergalactic dust is significant and could affect observations with NGST. 



Parameter 


Standard 


Range Considered 


Compton y c /10 5 i] y 


Opt. Depth T dust /r] T 




0.67 


0.67-1.0 


0.42-0.46 


0.40-0.69 


n 


1.0 


0.8-1.0 


0.38-0.42 


0.21-0.40 


IMF tilt (J3) 





-0.7-1.0 


0.17-0.54 


0.16-0.52 


H2 feedback 


yes 


yes /no 


0.42-0.46 


0.40-0.77 


n h 


0.05 


0.01-0.1 


0.02-1.71 


0.08-0.80 
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Fig. 1. — Top panel: the full spectrum at z=0 in our standard model for the star-formation 
history. The solid line shows the total spectrum, the dotted and dashed lines show the separate 
contributions of the original CMB, the stars, and the dust. The net contribution of dust results 
in absorption at A ^ 0.015 cm, and emission at longer wavelengths. Bottom panel: the redshift 
evolution of the fraction of baryons in stars, the dust temperature, and the y-parameter. Note that 
/star = -Fcoii^^^star/^b- The dashed curve shows the history of the CMB temperature. 
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Fig. 2. — Top panel: optical depth in our standard model to absorption and scattering by dust 
at z > 3, as a function of observed wavelength. The source is located at redshifts z s =10, 20, or 
100. Bottom panel: the expected flux with (lower curve) and without (upper curve) intervening 
dust, from star-clusters at z s =W and 20, containing 10 8 M and 4.6 x 10 8 M Q in stars. The stars 
are assumed to be distributed on the main sequence according to a Scalo IMF. The break in the 
spectra at the redshifted Lyman limit wavelength includes only the absorption by neutral hydrogen 
in the stellar atmospheres. 



